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We consider first-order phase transitions of the Universe in the radiation- 
dominated era. We argue that in general the velocity of interfaces is non-relativistic 
due to the interaction with the plasma and the release of latent heat. We study the 
general evolution of such slow phase transitions, which comprise essentially a short 



reheating stage and a longer phase equilibrium stage. We perform a completely an- 



alytical description of both stages. Some rough approximations are needed for the 
first stage, due to the non-trivial relations between the quantities that determine 



<]} ■ the variation of temperature with time. The second stage, instead, is considerably 

simplified by the fact that it develops at a constant temperature, close to the critical 
one. Indeed, in this case the equations can be solved exactly, including back-reaction 
on the expansion of the Universe. This treatment also applies to phase transitions 
mediated by impurities. We also investigate the relations between the different pa- 
rameters that govern the characteristics of the phase transition and its cosmological 



consequences, and discuss the dependence of these parameters with the particle con- 



tent of the theory. 
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1 I. INTRODUCTION 

a: 

It is well known that the Universe could have undergone several phase transitions in the 
early stages of its history, most of them associated with the spontaneous symmetry breaking 
of some symmetry. Some examples are the quark-hadron phase transition at the QCD scale, 
the phase transitions associated to the electroweak 577(2) x U(l) symmetry breaking or to 
Grand Unified Theories, and the Peccei-Quinn phase transition, related to the axion field and 
the strong CP problem. Cosmological phase transitions generically produce cosmic relics, 
such as topological defects, magnetic fields, or baryon number asymmetries, with potentially 
important cosmological consequences. The mechanisms for generating these relics build on 
the dynamics of the phase transition. 

In a first-order phase transition the dynamics is essentially determined by the nucleation 
and expansion of bubbles. At zero temperature, when a true vacuum bubble nucleates, it 
rapidly begins to expand with almost the velocity of light On the contrary, at high 
temperature, the bubble expands in a hot plasma, which is perturbed by the motion of the 
bubble walls. The plasma thus opposes a resistance to the expansion, that depends on the 
wall velocity. As a consequence, bubble walls feel a friction force, which prevents them to 
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accelerate indefinitely. Then, the velocity quickly reaches a stationary value, determined 
by the viscosity of the plasma and the pressure difference between the low temperature 
phase and the supercooled one. These quantities depend on the model, and it is known that 
the friction can be large enough to prevent the wall from acquiring relativistic velocities 
0, 0] Furthermore, the release of latent heat at the interfaces of the phase transition 
reheats the surrounding plasma up to a temperature that in most cases is close to the 
critical temperature. Consequently, the pressure difference that drives bubble expansion may 
decrease considerably, causing a drastic slowdown of the phase transition 0,0)- Therefore, 
at the radiation dominated era bubble walls generically undergo non-relativistic motion. 
Since the heat liberated at the interfaces is taken away by sound waves, the temperature 
can be assumed to be homogeneous, which simplifies the analysis. 

In order to have a first-order phase transition, the free energy must allow the coexistence 
of two phases. Therefore, we will assume that the free energy density T depends on some 
order parameter (usually a Higgs field) 0(x). In a certain range of temperatures, the free 
energy bears two minima separated by a barrier; one of them at = 0, which corresponds to 
the symmetric phase, and the other at a nonzero value m (T), corresponding to the broken 
symmetry phase. The difference in free energy density V (0, T) between some value of the 
order parameter and = 0, is generally given by a finite-temperature effective potential. 
The free energy difference between the two minima is thus given by 



At the critical temperature the two minima are degenerate, V (0 m (T c ) , T c ) = 0. Above 
the critical temperature the symmetric phase is the stable one, while below T c it becomes 
metastable, being m the absolute minimum. Finally, at some temperature To < T c , the 
barrier between the minima disappears, and the symmetric phase becomes unstable. At 
some stage between the temperatures T c and To, bubbles of broken-symmetry phase will 
be formed in the sea of symmetric phase. A bubble can be described as a configuration in 
which the order parameter is non- vanishing inside a spherical region (see e.g. |f|). After 
being nucleated, a bubble will grow with a velocity that depends on the pressure difference 
at the interface, Ap = —V (T), and on the viscosity of the hot plasma in which it expands. 

In this work we will be concerned with first-order phase transitions in the radiation- 
dominated epoch. Our aim is to study the development of such phase transitions within a 
completely analytical approach. Here we concentrate on the determination of the parameters 
that govern relic formation (e.g. number density of bubbles, bubble wall velocity, etc.). An 
analytic study of cosmological consequences will be addressed in JZ|, where the results of the 
present analysis will be used. In Section |H] we briefly review the influence of phase transition 
dynamics on the mechanisms for generating cosmic remnants in first-order phase transitions. 

In section IIHI we study a phase transition that completely develops at T = T c , with 
the two phases in equilibrium. This is a good approximation in the case of inhomogeneous 
nucleation in the presence of impurities. In the case of homogeneous nucleation of bubbles, 
this approximation correctly describes the evolution of the phase transition after some latent 
heat has been released. An interesting feature of phase equilibrium is that it is simple enough 
to solve analytically, including back-reaction on the expansion of the Universe. This is due 
to the fact that temperature is constant all the way through the phase transition. We thus 
can obtain the fraction of volume of the Universe that is occupied by the low-temperature 
phase as a function of time, with no need of any numerical calculations. 

Section IIVI is devoted to the analysis of the phase transition in the case of homogeneous 
nucleation. The main difference with the previous case is the initial stage of supercooling 
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and quick reheating back to the critical temperature. Section IS contains an analytical study 
of the phase transition dynamics. In sections IVII and I VI II we analyze the relations between 
the different physical parameters involved in the dynamics of the phase transition, leaving 
some technical discussions to the appendices. Our conclusions are summarized in section 

EM 



II. COSMOLOGICAL CONSEQUENCES OF PHASE TRANSITIONS 

Several cosmological objects may be formed in a phase transition of the Universe. Their 
abundance and characteristics depend on details of the development of the phase transition. 
Due to the complexity of the mechanisms by which these objects are created and the difficulty 
of describing the phase transition, several details of the dynamics (e.g. the variation of the 
nucleation rate or the wall velocity during the transition) are often disregarded for the sake 
of simplicity. An analytical study of the phase transition is thus important since analytical 
expressions will help taking into account the dynamics in a more rigorous way. In this 
section we review how phase transition dynamics affects the cosmological remnants. 

Electroweak baryogenesis. During the last two decades there has been much interest 
in the possibility that the electroweak phase transition could be the framework for the gen- 
eration of the baryon number asymmetry of the Universe (BAU). A first-order electroweak 
phase transition provides the three Sakharov's conditions for the generation of a BAU, al- 
though physics beyond the minimal Standard Model (SM) is mandatory in order to obtain 
a quantitatively satisfactory result (for reviews on electroweak baryogenesis see 0). Due 
to CP violating interactions of particles with the bubble walls, an asymmetry between left 
handed quarks and their antiparticles is generated in front of the walls of expanding bubbles. 
This asymmetry biases the baryon number violating sphaleron processes in the symmetric 
phase. The resulting baryon asymmetry is caught by the walls and enter the bubbles, where 
baryon number violating processes are turned off. 

It is important that the sphaleron processes be suppressed in the broken symmetry phase, 
in order to avoid the washout of the generated BAU when equilibrium is established after 
the phase transition. This requirement imposes a condition on the value of the Higgs field 
at the temperature of the transition cf) m (T t ) /T t > 1. Since 4> m is the order parameter, 
this is a condition on the strength of the first-order phase transition. 

The resulting BAU depends also on the bubble wall velocity. On one hand, if the velocity 
is too large, the left-handed density perturbation will pass so quickly through a given point 
in space that sphaleron processes will not have enough time to produce baryons; thus the 
resulting BAU will be small. On the other hand, for very small velocities thermal equilibrium 
will be restored and the baryon asymmetry will be erased by sphalerons; thus the BAU will 
be small again. As a consequence, the generated baryon number will have a peak at a given 
wall velocity, which is of order 10~ 2 |lfl[ 11, V^ . 



Both values of <p m and v w depend on T and vary during the phase transition. The 
dependence of the wall velocity is more critical, since reheating may cause it to descend two 
orders of magnitude before the transition completes. Baryogenesis may be either enhanced 
or suppressed by this effect jH [l3| , depending on which side of the peak the initial velocity 
lies. 



Baryon inhomogeneities. A general feature of cosmological phase transitions is the 
difference of particle masses between the high- and low-temperature phases. These mass 
differences give rise to different number densities in the two phases. At the QCD phase 
transition, for instance, baryons are much heavier in the hadron phase than in the deconfined 
quark phase. As the hadron phase expands, baryons are pushed into the quark phase 
region, leading to localized clumps of high density surrounded by large voids of low baryon 
density 

QEHI3- An important consequence is that large amplitude, small scale density 
fluctuations may survive until the nucleosynthesis epoch, affecting the standard scenario of 
big bang nucleosynthesis. Therefore, inhomogeneous nucleosynthesis may put constraints 
on the quark-hadron phase transition (see e.g. Moreover, if the quark phase reaches 

sufficiently high density, its pressure may balance that of the hadron phase. The quark 
matter trapped in small regions of space forms quark plasma objects that may survive until 
the present epoch [14 Il8|. 

Baryon inhomogeneities may also arise at the electroweak phase transition, since the 
amount of baryons produced through electroweak baryogenesis depends drastically on the 
wall velocity, and the latter has a considerable variation during the phase transition j^f. 
The geometry of the electroweak inhomogeneities is in general quite different from the QCD 
case. If the BAU peaks at a certain wall velocity, then the high density regions will form 
spherical walls, whose radius depends on the moment in the bubble evolution at which the 
peak velocity is attained. Furthermore, baryon number densities with the wron g sig n may 
arise in some regions of space, depending on the baryogenesis scenario 0, Eii llfljj. This 
gives rise to the interesting possibility of nucleosynthesis in the presence of antibaryons (see 
for example |2fj|). 

However, due to baryon diffusion and "neutrino inflation", baryon inhomogeneities gen- 
erated at the electroweak phase transition hardly survive until the nucleosynthesis time (see 
e.g. (MI1H)- Nevertheless, they may survive until the QCD scale. In that case, electroweak 
scale inhomogeneities can act as impurities for the quark-hadron phase transition (see next 
section). 

In summary, we can say that the amplitude and scale of baryon inhomogeneities generi- 
cally depend on the mean nucleation distance and on the variation of the velocity of bubble 
expansion. 



Topological defects and magnetic fields. If a global U (1) symmetry is sponta- 
neously broken at a first-order phase transition, the phase angle 9 of the Higgs field within 
each nucleated bubble is essentially constant, but phases in different bubbles are uncorre- 
lated. When bubbles collide, the discontinuity in the Higgs phase is smoothed out to become 
a continuous variation. The so called "geodesic rule" states that (for energetic reasons) the 
shortest path between the two phases is chosen (23j]. When three bubbles meet, a vortex (in 
two spatial dimensions) or a string (in 3d) may be trapped between them. This mechanism 
is obviously generalized to higher symmetry groups and other kinds of topological defects. 
Ignoring the dynamics of phase equilibration, it is easy to see that the number density of 
defects is proportional to the number density of bubbles. However, the final number of de- 
fects will depend strongly on the velocity of bubble expansion |24| . If the latter is much less 
than the velocity of light, then phase equilibration between two bubbles will have probably 
completed before they encounter a third bubble, thus reducing the chances of trapping a 
string. 

The above picture is in fact a rough simplification of the defect-formation problem. One 



complication is due to dissipation, since the Higgs field is coupled to the other fields in 
the thermal bath. Another complication arises when considering a gauge symmetry, and is 
caused by the fact that the phase of the Higgs field is not a gauge-invariant quantity, so 
it is convenient to define a gauge- invariant phase difference between two bubbles |24(. The 
phase difference is thus linked to the gauge field. In this case, dissipation can be taken into 
account by introducing the conductivity of the plasma. Then one can model the collision 
of three bubbles and calculate the evolution of the phase difference and gauge field. One 
can say that a vortex is formed whenever a quantum of magnetic flux is trapped in the 
unbroken-symmetry region between the three bubbles. 

From the above it is clear that the formation of local vortices is associated with the 
generation of magnetic fields. Therefore bubble collision constitutes also a mechanism for 
generating the cosmic magnetic fields (e.g., see |2fij). Of course, the magnetic field that is 
formed in this way corresponds to a spontaneously-broken symmetry which cannot be the 
electromagnetic U (l) em - Nevertheless, this mechanism can take place at the electroweak 
phase transition, where unstable cosmic strings 1 and hypermagnetic fields may be formed. 
The latter are subsequently converted to U (l) cm magnetic fields. It is interesting to note that 
the presence of magnetic fields may affect the dynamics of the electroweak phase transition 
(see e.g. 0). 

Calculating the magnitude of the magnetic fields and the density of defects that are left at 
the end of the phase transition involves the passage from three-bubble collision simulations to 
the computation of the full phase transition. Evidently, this is a difficult task. Although some 
simulations have been made (e.g., [281), several simplifications are generally required, which 
include forgetting about variations in the nucleation rate and the velocity of expansion of 
bubbles during the transition. An analytical investigation of phase transition dynamics may 
therefore clarify the picture and provide useful tools for the calculation of defect formation 
and magnetic field generation. 



III. PHASE EQUILIBRIUM 

We begin by considering the limiting case in which the first-order phase transition is as 
slowest as possible, namely, that of coexistence of the high- and low-temperature phases 
at the critical temperature TL Such a cosmic separation of phases has been studied for 
the QCD phase transition [l4 fl6l . liil |3(j . At the critical temperature there is no pressure 
difference between phases at the bubble walls, so the bubble expansion takes place almost 
in equilibrium. Assume at T = T c there are already regions with low-temperature phase. 
As the Universe expands, the fraction of space occupied by these regions increases, as the 
high-temperature phase converts to low-temperature phase. The loss of energy due to the 
expansion of the Universe is thus compensated by the latent heat released at the interfaces, 
and the temperature remains constant. In this scenario there is no supercooling. 

Since at T c the nucleation rate vanishes, such a first-order phase transition is only possible 
in the presence of impurities that induce the formation of bubbles. In this case inhomoge- 
neous nucleation theory applies. In a phase transition mediated by impurities there will 
still be some supercooling, which we neglect in this section for simplicity. The role of 
impurities in the early Universe could be played for instance by topological 0, or non- 



1 see for example |2j 



topological solitons 3^, 34], 35[. These may exist in the high-temperature phase, containing 
the low-temperature phase in their core. In this case their configurations become unstable 
or metastable below the critical temperature. When the system cools below T c these objects 
begin to expand and convert the Universe into the true vacuum. 

Another example of inhomogeneous nucleation is the case of the QCD phase transition 
in the presence of baryon number inhomogeneities ■ These may arise as a natural conse- 
quence of electroweak baryogenesis jfj, and can survive until the QCD scale |37f . Since the 
critical temperature is different in regions with different chemical potential [38|, |39| , bubbles 
will nucleate first in those regions with a higher T c i3| . If such regions are relatively small, 



and if they achieve the necessary amount of supercooling while the surrounding background 
reaches the critical temperature, then the baryon inhomogeneities operate as impurities for 
inhomogeneous nucleation. 

Even if the phase transition proceeds by homogeneous nucleation of bubbles, phase equi- 
librium will describe quite well a good part of the transition, whenever the latent heat is 
at least comparable to the energy density difference between the critical temperature and 
that at which nucleation effectively begins. In that case, the energy released will reheat the 
plasma back to a temperature very close to T c . At that moment bubble nucleation virtu- 
ally stops and the two phases remain close to equilibrium until the full latent heat of the 
transition is eliminated. We will analyze this case in the next section. 

The customary equation for the adiabatic expansion, p = —3H (p + p), tells us how the 
Universe takes energy from the hot plasma. Here p is the energy density, p is the pressure, 
and H is the expansion rate. It can be equivalently written in terms of entropy density, 
s = —3Hs, which is just the statement of entropy conservation, S = constant. It will be 
convenient to use it in the form s oc a~ 3 . At the beginning of the phase transition the whole 
Universe is in the symmetric phase, while at the end of the transition it is filled with the 
broken symmetry one, so we can write 

s = s u (T c ) a 3 /a 3 = s b (T c ) a 3 /a 3 , (2) 

where s u (p) is the entropy density of the unbroken (broken) symmetry phase, a is the cosmic 
scale factor, and = a (£»(/)) its value at the beginning (end) of the transition. We 
assume that, since the phase transition occurs very slowly, the latent heat released at the 
interfaces is quickly distributed throughout space and the temperature is homogeneous. We 
also assume that pressure and temperature remain constant during the phase transition. 
Phase coexistence at T = T c means that there are regions of space with different equations 
of state. Thus the entropy density has different constant values Sb (T c ) and s u (T c ) inside 
and outside the bubbles of broken symmetry phase respectively. The quantity s in Eq. (J2J) 
is the average entropy density of the whole system. The entropy in a comoving volume 
Vu = Vb + V u is the sum of two contributions, 5* = s b Vb + s u V u . Thus, 

s = s u + (s b - s u ) f b , (3) 

where fb is the fraction of space that is already in the broken-symmetry phase. The entropy, 
energy and pressure are derived from the free energy. At T = T c the pressure p c is the same 
in the two phases, Ap = —V (T c ) = 0. The latent heat of the phase transition is 

L = p u - Pb = T c (s u - s b ) = T C V (T c ) . (4) 

From Eqs. and (jH)), and using (J1J), we obtain the dependence of fb on the scale factor 



The dependence of the scale factor on time is given by the Friedman equation 



where for simplicity we have neglected the term k/a 2 . At constant temperature T c and 
pressure p c , the energy density is given by 

p = T c s-p c = T c s u a 3 /a 3 - p c , (7) 



where we have used Eq. (J2J) in the last equality. Inserting p in Eq. (JHJ and writing H 

3 



~a 3 da 3 /dt we can easily solve the equation for a 3 (t), 



where a> = ^QnGp c and the constant phase is determined by the initial condition a (tj) = a i; 
5 = arcsin ^p c /T c s u . 

During the phase transition we have two coexisting phases in the radiation dominated 
era, and we may consider different possibilities for the equations of state. The simplest one 
is to assume that the Universe is radiation-dominated before the phase transition, i.e., 

Pu = Pu/3, (9) 

with 

Pu (T) = tt 2 ^T 4 /30, (10) 

where g* is the number of effectively massless species 2 . In fact, this is not a realistic 
situation. In the symmetric phase the Higgs vev does not correspond to the true vacuum, 
so we should add a constant energy density to account for the energy of this state and have 
a negligible cosmological constant after the phase transition. It is interesting, however, to 
consider first this simpler case. Hence, at T = T c Eqs. (JHJ) and fTUj) give p c = p u /3 and 
T c s u = Ap u /3, so T c s u /p c = 4, 5 = n/6, and u = V3Hi/2. 

Before the phase transition H = ^87rGp u /3, with p u oc T 4 oc a -4 , so Eq. © gives the 
familiar relation H = (2t)^ 1 . The temperature descends like t _1//2 , and at T = T c the phase 
transition begins, since we are assuming that no supercooling occurs. Hence we can use the 
relation Hi = l/2tj as an initial condition. The scale factor thus takes the very simple form 




4 ^Mr-r - + 7T > ( n ) 



which does not depend on any parameter, i.e., during the transition the expansion of the 
Universe seems to be affected always in the same way, regardless the thermo dynamical pa- 
rameters of the theory. However, thermodynamics affects the duration of the phase transition 
and the subsequent evolution of the Universe, as we shall see immediately. 



2 In general depends on T. We will assume for simplicity that g* is constant during the phase transition. 
We discuss the effect of a variation of g* in Section IVII 



In the above equations it is apparent that the dynamics of the phase transition depends 
on thermodynamics only through the ratio 



r = 



T c s u 



(12) 



The fraction of volume is thus /& = r _1 [l — (a/a, 1 )" 3 ], with a/a, given by Eq. (JTTJ) . The 
phase transition concludes when f b = 1 [equivalently ]29|, when s in Eq. (J2J) equals sj, so 
its duration can be determined easily. Making use of some trigonometric algebra, 
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(13) 



Notice that Eq. ()13j) fails to give an answer for r > 3/4. The problem is that the scale 
factor given by Eq. ()11|) reaches a maximum when the argument of the sinus is tt/2. If r is 
small enough, that never occurs between £ f and tf. For r = 3/4 it occurs at t = tf, and for 
larger r it happens before tf, which means that the Universe begins to collapse before the 
phase transition has completed. This is not surprising. Indeed, if the energy density of the 
unbroken phase is given by Eq. (|10|). then the energy density of the broken symmetry phase 
is pb (T c ) = 7r 2 <7*T 4 /30 — L. So, there is a negative cosmological constant, which will begin 
to dominate sooner or later. If L < p u (i.e., r < 3/4), this won't happen during the phase 
transition, but below the critical temperature the Universe will collapse. 

In Fig. ^ we have plotted the expansion rate of the Universe from the beginning of the 
phase transition to the moment at which H becomes zero. The evolution of H without a 
phase transition is represented with a dashed line (long dashes). We have chosen a relatively 
large value r = 0.5 in order to get a visible departure during the phase transition, which 
occurs between t{ and tf (solid line). After the phase transition (short dashes) the expansion 
slows down, and the rate eventually vanishes at a time to- 
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FIG. 1: The expansion rate of the Universe for a phase transition at constant T = T c , in the case 
p u = 7r 2 <7*T /30 (negative cosmological constant). 

To consider a more realistic situation we must add a constant term ~ L to the initial 
energy density, so that we do not have a cosmological constant of the order of the scale of 
the phase transition. We take for simplicity 



p u (T) = tt 2 ^T 4 /30 + L. 



(14) 



It will be more convenient to re-express the general solution (JHJ) in terms of the conditions 
at t — tf, (a/cif) = (T c Sb/p c ) sin 2 [u> (t — tf) + 5'], with 5' = arcsin sJp c /T c Sb- From Eq (fT^ 
it follows that pb (T c ) = i^g^T^j '30, T c s u = 4/^/3, T c Sb 
Therefore we can write 



T c s u - L, and p c = p b /3 - L. 
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(16) 



where we have defined the time scale i = (2Hf) . Before the phase transition, p is given by 

Eq. (HID, so a oc sinh 1/2 ^- v /327rGL/3t). This has the form a ~ t 1 / 2 for t < (32ttGL/3) _1/2 , 

and departs from the radiation-domination behavior unless 7r 2 ^T 4 / '30 3> L. 
The fraction of volume in the broken symmetry phase is 
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and the duration of the phase transition is given by 
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It can be easily checked that, for small r, this solution coincides with the previous one. 
Furthermore, for r — > 0, the duration of the phase transition vanishes, as expected. At first 
sight, there seems to be a problem if L > pb/3 (i.e., for r > 1/4). However, all the previous 
expressions are still valid and real in the range 1/4 < r < 1. They can be written in the 
form 
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(21) 



For r — > 1 the duration of the phase transition becomes infinite because the constant energy 
density L is comparable to the energy density of the radiation, playing the role of a cosmo- 
logical constant that starts dominating at T ~ T c . Therefore the expansion of the Universe 
becomes too fast and the phase transition never ends. One expects that some of our initial 
assumptions will break down near this limit. For instance, the temperature and expansion 
rate will not be homogeneous, due to the significant energy density contrast between the 
two phases and the rapid expansion of the Universe. 



In Fig.|21we plot the expansion rate as a function of time for r = 0.8. We see that without 
the phase transition the Universe would enter exponential expansion (H — > constant). After 
the phase transition the evolution returns to the radiation-dominated relation H = 1/21 
By the end of the transition the departure of the expansion rate from its previous evolution 
becomes appreciable, because in this case tf — ti ~ t ~ ti. If the duration of the phase 
transition is short in comparison with the age of the Universe, the back-reaction on H can 
be disregarded. According to Eq. (fTH|) . this happens when the energy released is small in 
comparison to the energy density of the plasma (small r). 
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FIG. 2: The expansion rate of the Universe for a phase transition at constant T = T c , in the case 
Pu = TT 2 g,T 4 /30 + L. 



IV. SUPERCOOLING 

In the case of a phase transition mediated by homogeneous nucleation, bubbles start to 
nucleate at a temperature T < T c , when the gain in free energy inside a bubble is enough to 
compensate the cost of gradient energy at the surface. We have seen in the previous section 
that, even if bubbles begin to grow at T = T c , the phase transition may not come to an 
end if the parameter r is close to 1. The case of homogeneous nucleation is even worse due 
to the additional supercooling. As we shall see, a large latent heat is a general feature of 
strongly first-order phase transitions. In this case there may be extreme supercooling from 
which the Universe may never recover mug. Considerable supercooling and latent heat 
release may occur for instance in the quark-hadron phase transition ^|. Notice that in the 



case of large supercooling there will be an important departure from equilibrium, at least at 
the beginning of the phase transition. For the rest of the paper we will be mostly interested 
in the case r < 1. 

The nucleation and growth of bubbles in a first order phase transition has been extensively 
studied in the context of the QCD 0, and electroweak 0, II H H H H phase 
transitions. After a bubble is formed, it grows due to the pressure difference at its surface. 
There is a very short acceleration stage until the wall reaches a terminal velocity due to the 
friction of the plasma. It can be seen that this initial period in the history of the bubble 
expansion is negligible. We will assume again that the system remains close to equilibrium, 
in accordance with the assumption that the velocity of the bubble wall is small. If the wall 
velocity is less than the speed of sound in the relativistic plasma, v w < c s = \A/3, the wall 
propagates as a deflagration front. This means that a shock front precedes the wall, with 



a velocity v s h > c s . For v w <C c s , the latent heat is transmitted away from the wall and 
quickly distributed throughout space. We can take into account this effec t by considering a 
homogeneous reheating of the plasma during the expansion of bubblesja (For detailed 
treatments of hydrodynamics at different wall velocities see, e.g., |3. IsC ]). 



A. Bubble nucleation 

The thermal tunnelling probability for bubble nucleation per unit volume and time is 

mm 

r ~ A (T) e' Ss/T . (22) 

The prefactor involves a ratio of determinants associated with the quantum fluctuations 
around the instanton. In general it must be evaluated numerically. It is usually assumed 
to be roughly of order T 4 , since the nucleation rate is dominated by the exponential in 
(I22J) . We will consequently assume A (T) ~ T c 4 . S3 (T) is the three-dimensional instanton 
action, which coincides with the free energy of a critical bubble (i.e., a bubble in unstable 
equilibrium between expansion and contraction), 



S3 = Air / r dr 
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2 V dr 



+ V(0(r),T) 



(23) 



The configuration of the nucleated bubble may be obtained by extremizing this action. 
Hence it obeys the equation 

For temperatures very close to T c , the width of the bubble wall at the moment of formation 
is much smaller than its radius, and a thin wall approximation can be used 0, Hi| , in which 
S3 is expressed as a function of the critical bubble radius R c , the free energy difference V 
between the two minima of the potential, and the bubble wall surface energy a. The radius 
R c can thus be obtained by finding the maximum of S3. A similar approximation can be 
used to estimate the free energy and radius of a thick-walled bubble when temperature is not 
so close to T c |43j|. However, as pointed out in Ref. ji^l, due to the exponential dependence, 
the tunnelling probability may be strongly overestimated by using approximations to S3, 
leading to a sooner completion of the phase transition 3 . Since we do not intend to do 
numerical calculations in the present work, we will use the thin wall approximation. This 
approximation may be reasonable or not, depending on the amount of supercooling. 



B. Phase transition dynamics 

In the previous section we assumed constant temperature, and used Eqs. (j2J) and (JBJ) to 
obtain the fraction of volume fb in terms of the scale factor a. Then the Friedman equation 



In many cases the phase transition occurs in a tiny range of temperature about T c , so it is a good 
approximation to replace almost every quantity by its value at T ~ T c . The important exception are 
quantities such as S3, that depend directly on the free energy difference V (4>, T), which varies drastically 
with T at the critical temperature 0. 



determined a (t). In the present case the temperature is not constant, so we need an extra 
equation to solve for the three quantities /&, a and T. Such an equation arises by considering 
the nucleation and growth of bubbles 



h (*) = 1 - exp < — - 
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T(T')R(t',t) 3 dt ,[ j , (25) 

where T' is the temperature at t — t', t{ is the time at which the Universe reaches the critical 
temperature, 

U ^ iM P /Tl (26) 



where £ = a/907327t^^, and R (t', t) is the radius of a bubble that nucleated at time t' and 
expanded until the time t, 

RV,t)=B.m$ + f r ^$*r. (27) 

The factors of a in Eqs. (125)1 and ()27jl take into account the fact that the number density 
of nucleated bubbles is diluted, and the radius of a bubble enlarged, due to the expansion 
of the Universe from t' to t (see e.g. [47]). We can assume that this effect is negligible if 
the duration of the phase transition is small in comparison with the age of the Universe. As 
we have seen in the previous section, this is true when L/pb <C 1, which is the case we will 
consider. (This assumption will not hold, in principle, for the QCD phase transition). The 
wall velocity v w is determined by the equilibrium between the pressure difference V (T) and 
the friction force exerted by the plasma. The latter is proportional to the wall velocity. The 
constant of proportionality is the friction coefficient r\ (see section EH; so 

v w = -V (T) /n. (28) 

The exponent in Eq. (j2Hj) is minus the fraction of volume occupied by bubbles that nucleated 
between ti and t, if we do not take into account overlapping of bubbles. At the beginning 
of nucleation the formula /& (t) ~ (47r/3) J T (T") R (£', t) 3 dt' is correct. 
Using again Eqs. (0) and Q, we obtain the analogous of Eq. (J5J), 

h = m 1 7^ U (T) - s u (T c ) 4) , (29) 
s u (T) - s b (T) V a A J 

but since we already have an equation for f b , namely, Eq. (J23), we use Eq. (|29*Jl to express 
T in terms of fb and a, 

Eq. (13011 has come across within an approach that differs from previous works, so it is 
worthwhile spending a few words discussing its physical meaning. We may follow for instance 
Ref. 0|, and use energy (non-) conservation in the following way. On one hand, we may 
write the total energy in a volume Vu = V u + V& as E — [p u + Apfb] Vu = pVu, where 
Ap = pb — p u - If the Universe were not expanding, energy conservation during the phase 
transition would give 

p = p u + Apfb + Apfb = 0. (31) 



This gives the rate p u at which the plasma takes energy from the change of phase. On the 
other hand, when it is not undergoing a phase transition, the Universe takes energy from 
the plasma at a rate 

Pu = p = -Ap u H. (32) 
If we join the two equations, we obtain the total rate of change of energy as 

p u = -Apf b - Apf b - Ap u H, (33) 

from where we get an equation for T. However, if the phase transition and the expansion 
of the Universe are taken into account at the same time, additional terms appear both in 
Eqs. (|3*T|) and ()32|) . On one hand, there is a term of the form pVu in E, which produces a 
new term —3Hp. This just accounts for energy dilution. On the other hand, the expansion 
of the Universe takes energy not only from radiation, since it is not the only component in 
the equation of state. Bearing in mind the two coexisting phases during the phase transition, 
the rate at which energy conservation is violated is given by dE = —p u dV u — PbdVb- Using 
this to obtain p, we finally find 

p u = -Apf b - Apfb - 3H (p u + Pu + Apf b + Apfb) - A P f b , (34) 

where Ap = p b — p u . Since p + p = Ts, it can be seen that the first terms of this equation 
reproduce Eq. (j3HJ), but there are additional terms. Using s = dp/dT and rearranging 
Eq. (jnU), we see that it is just the equation for entropy conservation, 

s u + Asf b + As f b = -3H (s u + Ash) , (35) 

from which we may re-obtain the result (j3*Uj) . The discrepancies between Eq. f!33|) and 
Eq. (jMjl will not be important as long as the latent heat L is not significant and T remains 
close to T c . Indeed, we can neglect the last two terms inside the parenthesis in Eq. (jSlj), 
provided that Ap < p. The last term in Eq. (jMj) is responsible for the appearance of 
the entropy difference As = —V (T) as the factor of fb in Eq. ()35|) . instead of the energy 
difference Ap = V (T) — TV (T). Since V (T c ) = 0, the two quantities are related by 
Ap ~ TAs for T ~ T c . Therefore, Eq. (|33j) gives a good approximation in the case r«l. 
Still, the fact that in Eq. (J30|) the temperature is already integrated may constitute an 
advantage for analytical calculations. 

Finally, the evolution of the scale factor is given by the Friedman equation, 

H 2 = ^ [p u (T) + Ap/ b ] . (36) 

If Ap p, then one can use the customary formula H oc T 2 for the expansion rate. In 
fact, the variation of p u is of the same order of Ap, so if we neglect Ap, to be consistent 
we should also set H = const oc T c 2 in this approximation. Again, this is reasonable if the 
duration of the phase transition is short enough. 

In each particular model, Eqs. (j3Hj) and (j3Hj) can be used to calculate numerically the 
evolution of the phase transition. In this paper, dough, we will make analytical estimates of 
the parameters that are relevant for the cosmological consequences, such as the temperature 
and wall velocity in the different stages of the transition. 



C. Bubble coalescence 



When bubbles occupy more than 30% of space, they meet and percolate. This occurs 
when bubbles have a characteristic radius Rq ~ (0.3/n;,) 1 ^ 3 , where is the number density 
of bubbles. Bubble coalescence provides a different mechanism of bubble growth, in which 
the driving force is surface tension instead of pressure difference. When bubbles collide 
and percolate, they arrange themselves into a system of fewer, larger bubbles in order to 
minimize the surface area. We can estimate the characteristic time of this process as follows 
Assume that when two spherical bubbles of radius R meet, they form a single bubble 
of radius 2 1//3 i£. In this process, a mass of fluid m ~ pR 3 is moved a distance of order R. If 
this is done in a time t, the kinetic energy involved in the process is K ~ m (R/t) 2 . This 
energy is supplied by the surface of the bubbles. The surface energy released in the process 
is AE ~ aR 2 , so this occurs in a time t ~ {pR? / a) 1 ^ 2 . 

Once bubble coalescence begins, it will be the fastest mechanism of bubble growth if the 
rate t _1 is larger than the rate v w /R given by the wall velocity (|28jl. Therefore this process 
will dominate until bubbles reach a characteristic size R\ ~ a/v^p. If R± < Ro, then it 
never dominates. Otherwise, during the period in which the radius varies from Ro to Ri 
bubble coalescence is the fastest process and must be taken into account. 

The process of bubble coalescence may end in two different ways. If v w is large enough, 
the radius -Ri is very close to Ro, and after a short period bubbles continue to grow with 
velocity v w . If, on the contrary, v w is very small, then coalescence dominates for a larger 
period, until the low-temperature phase occupies more than 50% of the total volume. At 
this moment the regions of high-temperature phase detach into isolated bubbles and the 
process stops. The interface velocity is again determined by pressure difference, friction, 
and latent heat release, so these bubbles shrink with velocity v w . This occurs at a bubble 
radius R2 ~ (0.5 /rib) 1 ^ 3 . Notice that, although this expression is similar to that for Ro, in 
fact R2 may be very different from Rq since rib may decrease significantly during the process 
due to the expansion of the Universe. 

Finally, when bubbles of the high temperature phase are so small that surface energy 
becomes dominant over volume energy, the shrinking is accelerated until the symmetric- 
phase bubbles disappear or, eventually, until a topological defect is formed. 

V. BUBBLE NUMBER AND INTERFACE VELOCITY 

In the case of homogeneous nucleation, the Universe cools down to a temperature < T c 
before bubble nucleation becomes appreciable, so the phase transition effectively starts at 
a time > ti. Then, the phase transition proceeds in two main steps. At the beginning 
bubbles nucleate with a rate V (Tn), and expand with a velocity given by the friction coef- 
ficient 77 and the pressure difference V (Tjv). After a short time the energy released by the 
change of phase reheats the Universe up to a temperature T r close to T c . Then, a longer 
period begins, in which the two phases are close to equilibrium. How long is this period, 
and how close is T r to T c , depends on the latent heat. The free energy difference in this 
stage is V (T r ) ~ 0, so the velocity of bubble expansion decreases and the bubble nucleation 
rate becomes extremely suppressed. Therefore, this second stage may be very similar to the 
inhomogeneous nucleation phase transition discussed in Section ITTT1 

In order to avoid numerical calculations, we will need some approximations for the fraction 



of volume the temperature T, and the expansion rate of the Universe H [formulas (|25|). 
(jnOJ) and PS])]. We notice that the latter cannot be affected by the phase transition until 
the energy released becomes appreciable. Furthermore, when this occurs, the plasma has 
reheated up to T r ~ T c , and enters the phase equilibrium stage. We have seen in Section ITTT1 
that H is not modified significantly if the parameter r defined in f)12j) is small. So, we do 
not need to consider back-reaction on the scale factor and we will assume that the evolution 
of the Universe does not depart from the standard relations H ~ l/2t, a oc t 1 ! 2 . We will 
also need an approximation for the nucleation rate and free energy difference. 



A. Thin wall approximation and linearization of V 

If the width of the wall is much less than the radius of the bubble, we can neglect the 
second term in Eq. (}2*4*)) to obtain the wall profile (this approximation is exact at the critical 
temperature, at which R c — > oo, and gives the usual kink profile). If we then multiply 
by d(p/dr and integrate using the boundary conditions d<p/dr = and V = outside the 
bubble, we find that 

so, d(f)/dr = — \J2V since at the wall <fi falls from <p m to 0. Inserting this in Eq. (f2*3|) we 
obtain the free energy of the critical bubble in the thin wall approximation, 

S 3 = —RlV (T) + 47rR 2 c a (T) , (38) 

where the free energy difference V (T) is defined in Eq. (0) , and a (T) is the surface tension 
of the bubble wall, 

2 



° = j\Tr) dr = h V2Vd *- (39) 

Maximizing with respect to R c we get the values of the critical radius and action, 

R c = -2a /V, 
S 3 = 16na 3 /3V 2 . 

Since a does not change significantly during the phase transition, it can be approximated 
l-.v n (7; ). 

Since the thin wall approximation is valid when T does not depart significantly from T c , 
we can also make a linear approximation for the free energy difference, 

V(T)~L(T-T C )/T C . (40) 

So, the exponent in the nucleation rate (|2*2*j) becomes 

Ss(T) ^ W Tc 

T 3L 2 (T C -T) 2 ' V ; 

and the critical radius is 

R c ~2aT c /L(T c -T). (42) 
With the linear approximation (|4()|) . the wall velocity is given by 

v w = L (T c - T) /7]T C . (43) 



B. The first stage: nucleation and reheating 



Neglecting the effects of latent heat, the phase transition may be assumed to occur 
roughly at t = tjy, since it goes on for a very short interval St <C t^ ji^. With the inclusion 
of latent heat the phase transition evolves during a longer period, that begins at t^- In 
this case, we expect that the first part of the evolution, in which the Universe is reheated 
up to a temperature T r , has a time scale of the same order of the time interval St of the 
phase transition without latent heat. This is confirmed by numerical calculations ^3[. The 
nucleation rate T vanishes at T = T c , but it changes very quickly with temperature and 
becomes of order T 4 at T = T , where the barrier between the two minima of the free energy 
disappears. This is an extremely large rate, so it is impossible that the Universe supercools 
close to T [48]. We will assume that the temperature T N is close enough to T c that the 
approximations of the previous subsection can be used. 

In Ref. ji^] the onset of nucleation was assumed to occur when the probability that a 
bubble was nucleated inside each causal volume is one, 



The causal volume is given by Vh = d 3 H , where the horizon size dn scales like the age of the 
Universe, du ~ 2t. The cosmological scale t, however, is in general too large in comparison 
with the scale of phase transition dynamics, which at t = tjy is given by tjy — ti. The scale 
of phase transition dynamics is roughly determined by the temperature variation during the 
phase transition, which is bounded by the difference T c — Tq. 

Consequently, it may be more appropriate to consider a different causal distance d c , 
related to the dynamics of the phase transition in the following way. We may say that 
bubbles begin to "see" each other at a time tw when their mean separation is of the order 
of the distance travelled by a sound wave since time tj. Then, the causal distance is given 
by d c ~ c s (t N — ti), where c s = 1/V3 is the velocity of sound in the relativistic fluid. This 
defines a causal volume in terms of T/v 



We remark that the real improvement in using Eq. ()45|) instead of Vh does not come from 
considering the velocity of sound, which is ~ 1, but from the fact that in many cases the 
time elapsed from tj to t^ will be much less than the age of the Universe ti ~ ^Mp/T 6 2 . The 
volume V c is thus suppressed with respect to Vh by a factor [(Tj — T/v) /T] 3 . The nucleation 
time 1 7v calculated in this way is larger, since more bubbles need to be nucleated before they 
are separated by a distance at which they are causally connected to each other. 
There will be a bubble in each volume V c when 



To evaluate the integral in ()4fij) we make use of the following approximation [4j3] (see also 
30]). The three-dimensional action in Eq. (j41)l can be expanded about any temperature T, 
in the form 




(44) 



V c = [c s iM P {l/T 2 N -1/T 6 2 )] 3 , 



(45) 




(46) 




(47) 



where x = (T — T*) / (T c — T*). Since the integrand in Eq. (pftjj) is sharply peaked at T/v, we 
choose T* = T/v and use the expansion (|37jl to evaluate the integral. We find 

Q (rr \ irr /II 2 £Mp T C -T N S 3 (T N ) 

S 3 (T N ) /T N ~ 4 log — — + 4 log — log—; , (48) 

which will be in general dominated by the first term. For temperatures T several orders of 
magnitude below M P , we have S 3 (T N ) /T N > 100 (it is e.g. ~ 140 for the electroweak scale, 
and ~ 180 for the QCD scale). From Eq. (|41|) we obtain 

T c \ 2 3L 2 T C ( 2iM P 3L 2 T C T c - T N \ 

'41og— ^- + log- — ^ + 61og — , (49) 



T c — T N J 167rcr 3 \ T c 8%a 3 T c 

where we have used that T/v — T c . In some of the following estimations we will replace 
T c — Tjy by T c — T inside the log's, since both differences are generally of the same order. 

Immediately after t = tjy, the temperature increases at a rate which is given by Eq. 
(PH). Under the current approximations we can write p = Lf^ — ApH. Therefore, the rate 
of change of energy of the plasma is p/p ~ (3r/4)/ft — AH. Thus, on one hand energy is 
increased at a rate ~ rSfi/Sti, where Sfi is the fraction of volume converted to the low-T 
phase during the reheating stage, and 5t\ is the duration of this stage. The rate of energy 
decrease, on the other hand, is given by H ~ 1/1 Since 5t\ is much shorter than the age of 
the Universe, the increasing rate is much larger than the latter. The total change in energy 
density is thus 5p ~ Ldft. 

The temperature cannot increase beyond T c , so if L > p (T c ) — p (T/v), the fraction Sfi will 
be less than one. This means that before the phase transition completes, a final temperature 
T r very close to T c is reached, and the phase transition proceeds more slowly until 5f = 1. 
For L < p (T c ) — p (T/v) we have a variation 5p ~ L during the phase transition, which 
gives us an idea of the reheating that occurs. There may be significant reheating and a 
considerable variation of the wall velocity, but there will not be a long phase equilibrium 
stage, since in this case 5f\ ~ 1. The case L p (T c ) — p (T/v) corresponds to weakly 
first-order phase transitions, which most likely occur by spinodal decomposition rather than 
by nucleation and expansion of bubbles [H, 0]. We could also have L ^> p (T c ) — p (T/v), 
although we are assuming L p (T c ). In this case we have Sfi 1, which means that most 
of the phase transition happens with the two phases near equilibrium, and we can apply the 
analysis of Section IIII1 In this section we thus concentrate in the case in which the latent 
heat is comparable with the energy difference p (T c ) — p (T N ). 

The expansion of bubbles is governed by Eq. (|27)l . At T = T N bubbles nucleate with a 
radius R c given by Eq. (jl2j), and after a time 5t the radius has increased an amount v w St, 
with v w given by Eq. (|4*H|) . The ratio of the two distances is 

v w St L' 2 (T c -T N ) 2 St 

~W H 2noT- ■ (50) 

For the time scale of the phase transition dynamics, St ~ £Mp (T c — T ) /T 3 , we have 

v w St L\ M P fT c -T \ 3 

(51) 



R c rjaT T V T 



Thus, if the phase transition takes place at a temperature sufficiently below the plank scale, 
the bubbles will grow so rapidly that we can safely neglect the initial radius R c in Eq. (|2Tj) . 



At the beginning, bubbles expand with constant velocity v w (T N ). When reheating be- 
comes important, the bubble expansion slows down and we enter the second stage. Conse- 
quently, during the first stage f b can be roughly estimated without taking into account the 
liberation of latent heat. We can thus estimate the rate T/T ~ rf b without considering the 
back-reaction of the reheating on f b . It is only at the beginning and at the end of reheating 
that rfb ~ H and T ~ 0. Indeed, soon after t = tjy, the temperature takes its minimum 
value T m < T/v (See Fig. EJ), then it increases to T T < T c . The value of T m is important be- 
cause the nucleation rate turns on at t ~ tjy, is maximal at t = t m , and turns off again when 
the temperature has increased back to T/v. All the process occurs in a time 5tr ~ (t m — tjv), 
which is determined by the speed at which temperature changes. This interval Str in which 
T is not negligible is much less than 8t\, due to the exponential dependence of T on T c — T. 
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FIG. 3: Typical evolution of the temperature during a phase transition with supercooling. 



Accordingly, the temperature T m occurs when 



rAnv w (T m ) / T (T) v 2 w (T) {t m - tf dt ~ H. 



(52) 



where the time-temperature relation t = £M P /T 2 must be used to evaluate the integral. 
Using the expansion (}T7j) about T* = T m , we find 



S3 (T„ 
T 
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(53) 



The number density of bubbles is given by 



n b 



T (t) dt. 



(54) 



Since T (t) is sharply peaked at t m , we can estimate n b as 2 f* m T (t) dt and use again the 
dependence t oc T -2 . The integral can be calculated again using the expansion (|47j). which 
gives 



n b ~ T. 



, XM P T C -T„ 
Tc T c 



S 3 (Trr 

T 



n -1 



-S-i(T m ) /T m 



(55) 



We may define the interval Str by writing n b ~ V (t m ) Str- Then we can easily see in Eq. 
(1551) that 

St r r.[S 3 {T m )/T m )]- 1 (t m -t i ). (56) 

Therefore, t m —t^ is much less than t m —U (most typically about two orders of magnitude 4 ). 
Using the result (JHBJ), with T m ~ Tjy, we finally obtain the density of bubbles, 

"^^^{m;) (iwzvj T - (57) 

The bubbles thus nucleate during the short time Str about t = t m , and expand with a 
velocity V\ ~ f w (ijv) f° r a time until the temperature gets close to T c . According to 
Eqs. ipjl and <P|1 . 

1 /I67ra 3 \ 1/2 r ^ 1/2 

where K is a shorthand for the sum of logs in Notice that in Eq. (JSJ), v w oc L (T c — T), 
but also T c — T N oc L -1 , so f w only depends on L through the logs in if. 

If a non-negligible fraction of the volume is taken up by bubbles during the first stage, 
the interval Sti can be estimated from (47r/3) f 3 (t N ) St\n b ~ 1. This gives 



1/6 /g 3 (T ^" J/2 



c 



So, apart from a mo del- dependent factor, we find a general tendency to the relations Str "C 
<5ti < t/v — ^i- The value of 53 (T m ) /T m (and that of T m ) can be obtained similarly to the 
case of Tat- It is interesting to note that St\ has only a logarithmic dependence on the 
friction coefficient rj. This is because the dependence on the wall velocity is twofold. On one 
hand, the lower the wall velocity, the longer the time St\ needed to reheat the plasma. But 
on the other hand, the lower the wall velocity, the longer will also be the time Str m which 
bubbles are formed, and the larger their number. This causes a shorter St\, since there are 
more bubbles to produce the reheating. 



C. The second stage: phase equilibrium 

1. Inhomogeneous nucleation 

If the formation of bubbles is associated with the presence of impurities, the phase tran- 
sition occurs at T ~ T c , and the number density of bubbles n b is an external parameter that 
depends on the density of impurities. According to Eqs. ([15]) -111 7) 1 . for small r the evolution 
of the phase transition is given by 

f b ~ 3HT- 1 {t - U) , (60) 



4 This shows that it is more accurate to estimate the difference t m — ijv ~ 6tr in this way, rather than 
subtracting the values given by Eqs. gSJl and 



and the rate at which the phase transition goes on is fb = 3H/r, i.e., a factor of 1/r larger 
than the rate of expansion 3H of a comoving volume. Our assumption r < 1 implies that 
A > H and 5t < £. 

To calculate the velocity of the interfaces, we assume that all the bubbles begin to expand 
at t = ti. Thus, the fraction of volume occupied by bubbles is ft, = n b ^-Rf, where Rf, (t) 
is the bubble radius. At the midpoint of the transition we have fb ~ 4irnbR 2 v w , where 
R = (47rV3)~ 1/3 is the average radius. Therefore the mean velocity is given by v w = 

1/3 

(47rnb/3) H/r. We notice, however, that even in this case in which fb is constant, the 
wall velocity may change significantly during the transition, since v w oc R^ 2 . The total 
variation of v w depends on n&. 



2. Homogeneous nucleation 

In the case of a phase transition with supercooling, the situation is very similar after the 
plasma has reheated up to a temperature T r ~ T c . The transition proceeds at a rate 

f b ~47rv w (t) [ T(t')R 2 (t',t)dt'. (61) 

In this stage the nucleation rate has turned off. We have seen that T peaks sharply at a 
certain time t m in the previous stage, so we can write Eq. (JBTj) as fb ~ ^v^nbR 2 (t m , t), where 
t> 2 = v w (T r ), rib is given by Eq. (}5Tj) . and R (t m , t) ~ t> 2 (t — t m ). In any case, the integral in 
(JfiTj) is an average of the squared radius of the bubbles, and for the present estimations we 
can set ~ 47Tf 2 nb-R 2 , with R given by n b , just as in the case of inhomogeneous nucleation. 
Since the temperature is almost constant, Lfb — 4pH, i.e., all the released latent heat is 
taken away by the expansion of the Universe. Thus, again fb = 3H/r and the velocity 
coincides with that of the inhomogeneous nucleation case 

( 3 \ 1/3 H 

- (62) 

Although in Eq. (jo^j) it seems that the wall velocity during phase equilibrium does not 
depend on the friction, in fact it is proportional to r]~ l due to the dependence of rib- 
Since ~ H/r, the duration of this stage is 

8t 2 ~ rH~\ (63) 

The reheating temperature T r must be such that the pressure difference is adjusted so as to 
give the velocity (}62|) . Using Eq. (}43|) . we find 



T c - T r _ rjH_ ( 4nn b \ 
T c ~ rL V 3 J 



-1/3 

(64) 



As expected, the larger the latent heat, the closer will be T r to T c . Unlike 5t2, the values 
of T r and t> 2 are not easy to determine, since they depend on the number of bubbles that 
nucleated in the previous stage. Using Eq. (j37j) and taking into account that HM P /T 2 ~ 1, 
we find that roughly 

T c -T P ~f?^p^ (T C -T N ), (65) 



which confirms that generally, T c - T r < T c - T/v- According to Eq. (jjHj) , the same relation 
holds for the velocities v% and v 2 . 



In the range 0.3 < /& < 0.5 bubble percolation takes place. We have seen that this 
process gives a contribution to the bubble expansion rate, of order (a / pR 3 ) 1 ^ 2 . For /b ~ 1 
this rate is 



To establish the importance of this rate we should compare it with the rates StJ or St^ 1 , 
corresponding to the two stages we have studied. Such a comparison is difficult to carry out 
without specifying a model, so we will ignore this effect in the subsequent discussions. 

Although coalescence is bounded to occur in the above range of fb, it could have important 
consequences if the associated bubble growth rate is significantly larger than those given by 
the time scales 5ti and 5t 2 . In a specific model, comparison of Eqs. (j5H|) and with (ffiHj) 
should not be hard to do, once rib has been evaluated. 



We have seen that all the parameters that describe the dynamics of the phase transition 
(i.e. Sti, 5t 2 , etc.) depend on a few thermodynamical parameters, such as the latent heat 
or the friction coefficient. The formation of the different cosmological products of a phase 
transition thus depends on these quantities, and also on other parameters, such as the 
conductivity of the plasma. These quantities are physically related, since all of them come 
from the equilibrium or non-equilibrium thermodynamics of the same underlying theory. 
This should be taken into account in phase transition calculations, when ranges of parameters 
are considered. Unfortunately, it is hard in practice to establish general relations between 
these physical parameters. Depending on the theory, it may even be impossible to compute 
some of these quantities. 

One can gain some insight on the relations between thermodynamical quantities by con- 
veniently modelling the free energy. The problem is further simplified by referring to the 
general form of the perturbative effective potential. In that case the thermodynamical quan- 
tities can be related to the parameters of the microscopic theory. We dedicate this section 
and the following to study the aforementioned relations. We will concentrate only in those 
parameters which influence directly the dynamics of the phase transition. An analysis of 
other parameters that affect the generation of cosmological remnants is considered in Q. 



D. Coalescence 




(66) 



VI. THERMODYNAMICAL PARAMETERS 



A. Free energy and viscosity 



We assume the free energy density takes the form 




(67) 



where the scalar field is the order parameter, and 




(68) 



is the free energy density difference between the symmetric and the broken-symmetry phases. 
The parameter is the number of light species of the plasma. In general, g* depends on 
temperature, but it is usually approximated by 

7 . 

bosons fermions 

where the sums are on particles with masses < T, and gi is the number of degrees of 
freedom of species % (see appendix IXJ) . 

It depends on each particular case wether Eqs. (fHTj) and can be derived from the 
microscopic theory. At any rate, they can be regarded just as a simple model for studying 
the dynamics of the phase transition, being the latter first-order if the coefficient E is 
nonvanishing. The parameters of V ((f), T) can be chosen in such a way that the free energy 
carries the thermodynamical properties of the theory we wish to study 0, HI 0]. For 
instance, these parameters determine the values of the critical temperature, latent heat, 
surface tension, and correlation length. The thermodynamical parameters could be obtained, 
e.g., with lattice simulations (see for example [50]). Then one can use those values to 
calculate the parameters To, D, E and A. Furthermore, in general the order parameter is 
a Higgs field or a combination of Higgs fields, and V (0, T) is the finite-temperature effective 
potential (see e.g. We will consider this case in the next section. 

The effect of viscosity on the propagation of the bubble wall is calculated by considering 
its equation of motion in the hot plasma, 

□, + m+ £a*£/_^_ /l(t , l)=0 , (70) 

which can be derived by energy conservation considerations jH 0, 0, 0] • Here V (0) is 
the zero temperature effective potential, the sum is over all particles that couple to 0, m; 
are the 0-dependent masses (see the appendices), and are the phase space population 
densities. This equation can be obtained by thermally averaging the operator equation for 
0. If we separate / into the equilibrium population f plus a small deviation 5f, we obtain 
the equation 

where V (<j>, T) is the finite temperature effective potential, given by Eq. (|B5|). Since the 
departure from equilibrium is proportional to the velocity of the bubble wall, it is the last 
term in (ffT|) which gives the friction force of the plasma. 

A simple approach to the calculation of the wall velocity 0, H| consists in replacing the 
last term in Eq. (j71|) with a typical damping term of the form d<p/dt. Due to Lorentz 
invariance this term must be in fact of the form u^d^tfi, where is the four-velocity of the 
plasma. Eq. (fTTj) then may be written as 

□0 + V + (rjT) u^d^ = 0, (72) 

where fj is a dimensionless damping coefficient that depends on the viscosity of the medium. 
Boosting to a frame that moves with the wall, and assuming stationary and non-relativistic 
motion in the z direction, we have 



<P" = V (0)-^0', 



(73) 



where 0' = d(j)/dz. Multiplying both sides by 0' and integrating over — oo < x < oo we 
obtain 

fjTav w = V(T), (74) 

where V (T) is the free energy difference between the two phases, defined in Eq. (0), and a is 
the surface tension of the wall, given by Eq. (|39[). We have assumed here that temperature 
is constant across the wall. This is right if the wall velocity is small enough, so that the 
latent heat it releases has time to be uniformly distributed throughout space. 

Hence, the pressure difference is equilibrated by a friction force proportional to the wall 
velocity. The constant of proportionality is the friction coefficient rj = fjTa. Since the 
tension of the wall is related to the wall width L w by o ~ 0J^/ L w [see the first integral in 
Eq. (JHHJ)], we can also write r\ = fjTcj)^/ L w . 

A shortage of modelling the viscosity of the plasma in this way is that f) is a free parameter. 
The correct expression for r\ can be derived from Eq. (|7ip. In appendix El we show that 
particles with a thermal distribution give a friction coefficient 

Vth ~ fj th ^a, (75) 

while the contribution of infrared gauge bosons is 

T 3 

Vir — Vir-j—- (76) 

Evidently, both formulas agree with the above result if m ~ T. This treatment also allows 
for the evaluation of the coefficients fj, which depend only on the particle content of the 
plasma. 



B. Thermodynamical quantities and phase transition dynamics 



The free energy given by Eqs. ()67j) and bears a first-order phase transition, with 
two minima separated by a barrier. The critical temperature is related to T by 



rp2 rp2 

J^2 



\D' 



(77) 



At T > T c the global minimum of the potential is (f) = 0. At the critical temperature the 
two minima become degenerate, and below this temperature the stable minimum is 



0m (T) 



3ET 
~2\ 




(78) 



At T = T the barrier between minima disappears and = becomes a maximum of the 
potential. Therefore the phase transition occurs at some stage in between T c and T . The 
value = corresponds to the symmetric, high-temperature phase, and ^ corresponds 
to the broken-symmetry, low-temperature phase. The jump of the order parameter from the 
high temperature phase to the low temperature one is thus 



m (T c ) = 2ETJX. 



(79) 



According to Eqs. (|57jl and (jSHJ), the free energy density of the symmetric phase is 

fu = -^g*T A + L. (80) 

This gives the equation of state of a hot relativistic plasma with a positive cosmological 
constant 

Pu = —g*T A + L, Pu = Pu /3 - 4L/3. (81) 

The free energy of the broken-symmetry phase is T = T u + V (T) . The energy density of 
the broken-symmetry phase is pb = p u + Ap, with 

Ap = V (T) - TV' (T) . (82) 

The entropy density of the symmetric phase is s u = 2ir 2 g*T 4: /45, and that of the broken- 
symmetry phase is Sf, = s u — V (T). The latent heat of the phase transition is given by 
L = Ap (T c ) = T c As (T c ); hence, 

L = 8D^ T 2 C T 2 Q . (83) 

Comparing with Eq. (JZHJ), we find the relation L = 1D<$^Tq between the discontinuity of 
the order parameter and that of the energy density. As expected, strongly first-order phase 
transitions (i.e., with large <p m ) have large latent heat. 

The surface tension of the bubble wall is given by Eq. (jHSj) in the thin wall approximation. 
At the critical temperature the effective potential is given by 

V(<P,T c ) = ^^x 2 (l-x) 2 , (84) 
where x = \<p/2ET c = 4>/4> m . Hence, Eq. (J3HJ) is easily integrated and 

Although we haven't used it explicitly, in this approximation the field configuration near 
the wall can be solved analytically with the help of Eq. (|84jl . and gives the kink profile, 



^W = y 1 +tanh— ), (86) 



where 

L w = tfjZa (87) 

is the wall width 5 . 



5 L w may change during the bubble expansion due to the friction with the plasma |2| • We shall neglect this 
effect. 



Using Eqs. (|79|). (|H3)l and (|87j) we find the values of the friction coefficients (175} and (|75j). 



8y^E 5 r/ th 4 

% = ^T c \ (88) 



The two contributions have different parametrical dependence, so each will dominate in 
different regions of parameter space. For instance, if E <C A the infrared boson contribution 
may be much larger than that of thermal particles. The maximum velocity of bubble walls 
occurs at T ~ TV- According to Eqs. (JH%|) and (JHEJ), this velocity is the smallest among 

\3/4 p7/2 

~ ^-^" 1/2 > ~ tStI-^" 1/2 - ( 89 ) 
E 1 /2 ?7th A 13 / 4 ?7 ir 

To determine which one is the correct, it is necessary to know the relations between the 
coefficients E, A, and fj. We see that in the opposite limiting cases E <C A and E ^> A, one 
of the two velocities is <C 1, unless fj is too small. In the case E ~ A, the wall velocity will 
be small if one of the conditions E 1 ^ <C 77th or E 1 ^ <C fj- n is fulfilled, which does not seem 
unlikely in general (see the next section). This supports the assumption of non-relativistic 
wall velocities. 

It is evident that with the aid of the model ()67|68jl we can get more information about the 
generalities of phase transition dynamics. For instance, if we write Eq. (|49|) as a function 
of E, D, and A, and compare with Eq. (j77|). then we can locate the nucleation temperature 
in the interval T c — Tq, 

T C -T N ~ fvT^ 1/2 {Tc - T ) , (90) 

If E and A are comparable, this gives a value of T c — T/v roughly an order of magnitude less 
than T c — T . 

The relations between the different quantities that determine the dynamics of the phase 
transition are apparent in the above expressions. Specific relations will be of interest for 
different cosmological consequences. As an example, let us consider the effect of modifying 
the theory in order to obtain a more strongly first-order phase transition. To do that, we 
have to enlarge the value of the order parameter. Assume we accomplish this by increasing 
the value of the parameter E and keeping the other parameters invariant [see Eq. (j79)) ]. 
Then, there will be more supercooling, and one expects a larger departure from thermal 
equilibrium, since the pressure difference at T = T/v will be larger. However, according to 
Eqs. (|83p and (|88p. L and rj also increase. This tends to decrease the wall velocity in the 
two stages of the phase transition, in opposition to the effect of supercooling. 

In this work we assume for simplicity that g* remains constant throughout the phase 
transition. In fact, the number of effectively massless degrees of freedom may change during 
the phase transition. It is conceivable that some particles acquire large masses and decouple 
from the thermal bath; then Ag* = g m — g^ > 0. For instance, during the quark-hadron 
phase transition g* changes substantially. It is interesting to note that such a change may 
affect considerably the dynamics of the phase transition, even in the case A g* <ti g*. The 
effect of a decrease of g* during the phase transition is twofold. To begin with, the free 
energy of the broken-symmetry phase is larger than in the case of constant g*, so the critical 
temperature is lower (it is given by V (T c ) = — n 2 A g^T^/ 90 [13,|4l|). Therefore the phase 



transition is stronger 6 , and the latent heat T C V (T c ) is larger. In addition, the entropy 
released by the decoupling species gives an extra contribution of 4n 2 Ag*T c 4 / '90 to the latent 
heat. This contribution is comparable to the value of L as given by Eq. (|83p. if Ag* > 
D (E/\y (1 - E 2 /XD). In the case of a perturbative effective potential, this condition may 
be easily fulfilled for A ~ 1. 



VII. THE PHYSICAL QUANTITIES IN PERTURBATION THEORY 

If perturbation theory is applicable, the one-loop effective potential at high temperature 
has oftentimes the form of Eq. (jfi8j) . with parameters generally given by 

__ \ ^ 9ih 2 \ ^ gjhf 
24 ^ 48 ' 

bosons fermions 

DA' [J } 

3 ^ 127r' 

gauge 
bosons 

A = m 2 h /2v 2 , 

Here, hi are the couplings of the particles with <p, rrih is the Higgs mass, and v its zero 
temperature vev. The coefficient E in general involves only gauge bosons. In appendix 
1X1 we review the derivation of these results and discuss on the general assumptions and 
approximations that lead to Eqs. (}68|) . ()69)1 . and (}9T|) . In the discussions that follow we 
will sometimes take the electroweak theory as a reference point. The parameter To gives the 
temperature scale of the phase transition. Its order of magnitude is determined by m^, so 
it may be quite less than the scale v if A is small. Anyway, for the dynamics of the phase 
transition, the difference T c — T is more important than the temperature scale T . 

Regarding the viscosity of the plasma, we show in appendix El that the contribution of 
thermal particles to the parameter fj is given by 



Vth 



2 



where Xi — 2 for fermions and Xi — h% f° r bosons. Therefore the contributions of bosons to 
7] have an enhancement of (log h^ 1 / log 2) 2 with respect to fermions with the same Yukawa 
coupling. For instance, for h ~ 0.1 the boson enhancement is ~ 10. This means that friction 
may be much stronger in supersymmetric theories than in non-supersymmetric ones. For 
instance, it was found in Ref. [3J that a light stop may slow down the electroweak bubble 
wall in the MSSM an order of magnitude with respect to the SM. The enhancement is larger 
for lighter particles, but these do not contribute to the friction due to the h\ dependence. 
The contribution of infrared gauge bosons is 

fju ^ V log [m b {(j) m ) L w \ . (93) 



6 This could be important for baryogenesis 



Here, the sum is only on gauge bosons, but the coefficient g also involves a sum over particle 
species (see appendix |Bj) . Furthermore, the gauge coupling appears only squared, which 
means less suppression. The log enhancement in this case is ~ log (hb^/a) ~ (log /if,A -1 / 2 ) . 

It is important to compare the value of the parameter E with the other parameters, 
since E is responsible for the first-order nature of the phase transition. We can see in the 
formulas of the previous section that all the thermodynamical quantities are proportional 
to some power of E, while the parameters D and A usually appear in the denominators. In 
the perturbative approach ([91)1 . this parameter is generally smaller than the others. This 
is because E is a sum of gauge couplings to the third power, weighted with gauge boson 
degrees of freedom, while D involves squared couplings, and the sum is over all degrees of 
freedom. Regarding A, it can be comparable to E, but this constrains the value of the Higgs 
mass. 

The smallness of E indicates a tendency of perturbative effective potentials to give weakly 
first-order phase transitions 7 . This is apparent in the dependence of the order parameter, 
4> m /T ~ E/X, or in the temperature interval in which the first-order phase transition can 
occur, (T c — To) /T c ~ E 2 /XD. For example, in the case of the electroweak phase transition, 
we have E ~ 10~ 3 and D ~ 10 _1 for the minimal Standard Model. If we take a non- 
realistic value for A ~ E to get an order parameter of order T, we find a temperature 
range T c — T ~ 10 _2 T C . In the specific case of the electroweak phase transition, a small 
value of the Higgs field m (T) is undesirable for electroweak baryogenesis. In general, if 
this parameter is too small, the perturbative approximation breaks down (see appendix IA"|) . 
In the electroweak case, the way out is to consider extensions of the SM which provide 
additional bosons that contribute to the parameter E (HH, 0| • 

In the previous section we found that if we increase the parameter E while keeping the 
others constant, then we get a stronger phase transition and larger supercooling, but also 
larger values of 77 and L, which slow down the dynamics. It is evident that if E is augmented 
by adding particles to the theory, the value of D enlarges too, giving an additional increase 
of the latent heat L. If we add only a boson, the relative change will be more appreciable 
in E, because there are only a few terms in its expression, but if the boson comes together 
with several new species (as in the case of supersymmetry) , then the change of D will be 
much more substantial. According to Eqs. (jQZ|) and (jHSJ), the friction coefficient will also 
increase significantly when adding bosons to the theory. 

If E 2 /\D < 1, then T c ~ T , and L/T c 4 ~ 8D (E/X) 2 . It is interesting to compare the 
value of L with that of 5p = p (T c ) — p (T N ), to assess the effect of reheating, as discussed in 
section El Using Eqs. (j9*Uj) and (|77j). we may write 

SP - ~ K ~ 1/2 ^I- ( 94 ) 



Therefore, 

L ( 30K 1 / 2 \ / D 2 



L ( 30KV 2 \ f D 2 \ 



(95) 



The first factor is likely of order 1 and depends essentially on the energy scale of the transi- 
tion. The second factor is determined by the dynamics. It depends mainly on D, and may 
vary considerably if we change the particle content of the theory. Exemplifying again with 



7 We are only discussing one-loop order here. Things may be different at two loops |4: 



the electroweak theory, D can vary from ~ 10 _1 in the SM to D > 1 in the MSSM, so we 
pass from little reheating in the first case to large reheating in the latter. We remark that 
things may be quite different if Eqs. (|91|) are not valid. For instance, in the case of the 
quark-hadron phase transition L and 5p are typically of the order of the energy density p. 



VIII. CONCLUSIONS 

In this paper we have performed an entirely analytical study of first-order phase tran- 
sitions in the radiation-dominated era. We have seen that typically the high-temperature 
phase is supercooled to a temperature T/v, after which the transition proceeds in two steps, 
as sketched in Fig. El The first stage is complex, and some rough approximations must be 
made for an analytical treatment. Nevertheless, it can be checked with numerical results, 
(e.g. that the orders of magnitude are correct. The second stage is much more simple, 

since bubble nucleation has effectively stopped and bubbles expand very slowly. This stage 
develops very close to the critical temperature, with almost zero pressure difference between 
the two phases. Therefore, this part of the evolution is similar to the case of inhomogeneous 
nucleation, in which the presence of impurities induces bubble nucleations without need of 
supercooling. 

We have studied the case of a phase transition at phase equilibrium in some detail, taking 
advantage of the fact that it can be solved analytically for any value of the relevant parameter 
r = L/T c s(T c ) ~ L/p(T c ). This approach has allowed us to calculate the back-reaction on 
the expansion rate H, which is important for large r. It is well known that supercooling may 



lead to exponential expansion of the Universe [32, |40| . This possibility has been considered 
not only in the context of inflationary models, but also for the quark-hadron phase transition 
j47t | . Although our approximations break down for L > p(T c ), we observe the manifestation 
of the energy of the false vacuum for large L. Even if the phase transition begins at T = T c , 
when L is comparable to the energy density of the plasma the transition may take a long 
time to complete due to vacuum energy dominance. 

For the more probable case of a phase transition with variation of temperature, we have 
given a derivation of the integro-differential equations that govern the dynamics. In par- 
ticular, we have found a simple algebraic relation between the temperature, the fraction of 
volume occupied by the low-temperature phase, and the scale factor of the Universe, which 
holds under the usual assumption of adiabatic expansion. Using the thin wall approximation 
and the linearized form of the free energy difference, we have found analytical formulas for 
all the quantities that characterize the dynamics of the phase transition and may be relevant 
for the determination of its cosmological consequences. These parameters are the durations 
5ti and St2 of the two stages of the transition, the wall velocities vi and V2 at each stage, the 
total number density of bubbles, the time interval Str in which bubble nucleation is active, 
etc. We have expressed these quantities in terms of those that determine the dynamics, 
namely, thermodynamical parameters like the latent heat L, the wall tension a, or the fric- 
tion coefficient 77. As expected, for the phase equilibrium stage we have simple expressions, 
of the sort St2 — rH^ 1 , with obvious interpretation. More complex formulas arise instead 
for the reheating stage. Although these parameters must be calculated in each particular 
case, some relations can be established, that hold quite broadly. They allowed us to confirm 
some natural premises for the dynamics, e.g., that 5tr <C 5t\ <C <5t 2 and V2 <C t>i. 

We have studied also the interrelations between the thermodynamical parameters. When 
necessary, we made use of a simple model for the free energy. Aside from reproducing the 



desired features of the phase transition, it is well known that this model corresponds to the 
simplest high-temperature effective potential that arises in perturbation theory. We have 
derived general expressions for the parameters of this potential, which is useful in estab- 
lishing further relations between the thermodynamical parameters. We have also derived 
general expressions for the friction on the bubble walls. This is caused by the perturbation 
from equilibrium of the particles of the plasma due to the motion of the interfaces. We have 
compared the case of thermal particles to that of coherent infrared bosons. We have found a 
different parametric dependence of each contribution, which indicates that each of them will 
dominate in different parameter ranges. We have argued that probably one of these contri- 
butions will cause the wall to move non-relativistically. This justifies the near-equilibrium 
approximations that simplify the analysis of the phase transition. 

We have thus been able to find some general relations between the parameters that 
determine the dynamics. For instance, we have seen that if the first-order phase transition 
is strengthened, then the supercooling is intensified, but also the friction and latent heat are 
generically enlarged, giving a slower evolution. This general feature is easily detected with 
ad hoc variations of the parameters of the free energy, and further confirmed by the relations 
that arise for perturbation-theory values of these parameters. The amount of supercooling 
is characterized by the difference between the energy density of the plasma at the critical 
temperature, and that at which nucleation begins, 5p = p (T c ) — p (T N ). The specific relation 
between dp and the latent heat is decisive for the phase transition dynamics. As we have 
seen, the ratio L/5p can either be small or large, depending on the theory. Of course, it is 
larger for stronger phase transitions. We have argued that the case of interest corresponds to 
a latent heat comparable to S p. On one hand, a small L is related to too weakly first-order 
phase transitions. On the other hand, for large L most of the phase transition occurs close 
to phase equilibrium, and can be described as a phase transition at constant T ~ T c . This 
includes the case r ~ 1 if 5p <C p, thus justifying the approximation r < 1 in the case of 
supercooling. 

For the study of the different cosmological consequences, additional specific relations will 
be relevant in each case [?} , which can be obtained from the present analysis. Our analytical 
approximations will thus prove useful to include details of phase transition dynamics in the 
calculation of cosmic remnants, particularly with regard to the variation of the pertinent 
parameters throughout the phase transition. For example, the importance of the phase 
equilibrium stage has been already investigated in (H, in the context of electroweak 
baryogenesis. Other cosmological consequences are affected as well. For instance, the fact 
that the nucleation rate is turned off due to reheating evidently modifies the number of 
nucleated bubbles, and thus the density of topological defects. 
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APPENDIX A: PERTURB ATIVE FREE ENERGY 



Following |6(, we will obtain the high-temperature effective potential (or free energy) 
in the one-loop approximation, including leading-order plasma effects. Additional terms 
appear at higher-loop order. For instance, potentially important terms of the form 2 log 
arise at two loops |42| . However, inclusion of two- loop corrections makes the situation more 
complicated and lies out of the scope of our general analysis. 

We consider a gauge theory which is spontaneously broken by a vev of a scalar field 0. 
The tree-level potential is 

W) = -^V + ^ 4 , (Al) 

so the vev is given by v 2 = m 2 /\ and the Higgs mass by m\ = 2Xv 2 . 

At one-loop the effective potential picks up zero-temperature and finite-temperature 
corrections. With a cut-off regularization and tree level values for v and m^, the zero- 
temperature contribution of a particle species is |43[ 



' m 2 (0) 
log 



m 2 (v) 2 



2m 2 (v) m 2 (0) \ , (A2) 



where the ± is for bosons (fermions), g is the number of degrees of freedom of the species, 
and m (0) is the mass of the particle in the presence of the background scalar field. It is in 
general of the form 

m 2 (0) = fi 2 + h 2 (f) 2 , (A3) 

where h is the coupling of the particle with (i.e., the Yukawa coupling, gauge coupling, 
etc.). The finite-temperature corrections are of the form 

_L 9 J / 2 



"2tt 2 ^ 



J B , F (m 2 (0)/? 2 ), (A4) 



where the functions Jb and Jp can be expanded in powers of m/T for m < T, and fall 
off exponentially for large m/T. Therefore species with m ^> T decouple from the plasma 
and we make a high-temperature approximation in which we consider only particles with 
m < T. Expanding up to O (m/T) 4 we have (see e.g. p) 



^ ^4 m h rl 



T (0, T) = constant - ^9*T A - -^(f) 2 (A5) 

9b m b (0) 



127T , 

O 



+ V - T A log ( ^t) mi (0) + Y log 

b \ b/ f 



T 2 A f 



m f (0) 



where g* is the effective number of relativistic degrees of freedom, given in Eq. ()69j) and 
rrii = mi (0 = v ) are the physical masses. We assume that particles contributing to JF do 
not decouple during the phase transition, i.e., that the condition m (0) <C T is preserved in 



the range of temperatures of interest. This is a reasonable assumption provided that <C v 
for temperatures close to T c , which is consistent with the one-loop approximation. If some 
particles decouple from the plasma during the transition, the main effect is a change in g>*. 

The m 3 term is the contribution of the bosonic zero modes to the one-loop effective 
potential. This term is the most important to us; without it the phase transition would be 
of second order. However, for the zero modes the loop expansion has an infrared problem. 
The perturbative expansion breaks down at higher-loop order, since higher loops contribute 
powers of a = h 2 T 2 /m 2 and of /3 — h 2 T/m((j)) [6|. The way out is to dress the zero 
modes with daisy and superdaisy diagrams. The result of this resummation is a contribution 
of the form (to all order in a and to O (/?)) 

b 

plus contributions proportional to 2 which are unimportant within the present approxima- 
tions. Thus the mass gets replaced with its Debye mass 

Ai 2 = m 2 (0) + n b (0,T), (A7) 

where 11^ is the self-energy of the boson particle. In general it is a combination of squared 
coupling constants times T 2 . The exception are the transverse components of the gauge 
bosons, for which II = 0. 

If we replace the Debye mass (|A7|) in the cubic term of Eq. (|A5J) . and the masses (|A3j) 
everywhere else, the resulting terms can be grouped as follows: 

a. Constant terms. These are contributions to the cosmo logical constant. The total 
cosmological constant must be set by hand, so that it is almost zero after the phase transition. 

b. T -dependent, <fi -independent terms. Apart from the first term in (|A5[) . there are also 
T 2 and logarithmic terms, but these are of order (/ij/T) 2 and (/ij/T) 4 with respect to the 
T 4 term, so we can neglect them within the approximation <C T. We notice, however, 
that we could have a large negative //j of order T, such that r/ij is small (see below). In any 
case, these corrections modify the equations of state of both phases in the same way, and 
we do not expect them to affect significantly the dynamics of the phase transition. 

c. cf) 2 terms. The coefficient of eft 2 is the sum of a term proportional to T 2 , a constant 
term m\/A, and other constant and logarithmic terms, which are ~ h 2 m 2 /32n 2 . The latter 
are suppressed unless the Higgs mass is too small, so we will disregard them. In any case, 
these corrections are inconsequential for our purposes. They contribute to the value of the 
characteristic temperature T of the phase transition. However, for the dynamics of the phase 
transition the precise value of T is not relevant; the important parameter is the relative 
difference between this temperature and the critical one, (T c — T ) /T , which is independent 
of T . 

d. 3 terms. The A4 3 term has contributions from all the bosons, proportional to 

(Z^ + ^ + n,) 372 . (A8) 

These terms may strongly affect the nature of the phase transition, depending on the value 
of + lib (T c ). There are two limiting cases, 

i) if /i 2 + (T c ) ~ 0, (jA8|) contributes a term of the form T<f) 3 to the free energy, which 
favors a strongly first-order phase transition, and 



ii) if + II;, (T c ) ^> hl<p 2 , we can expand (|A8|) in powers of 0. This gives higher-order 
corrections to the coefficients of (p 2 and 4> 4 , and no contribution to 3 . 

Almost all particles fall in the second case, since II is of order h 2 T 2 , and ^ C T near 
the phase transition. The transverse components of the vector bosons, on the contrary, are 
protected against this thermal screening. Another exception may be a scalar with a negative 
\i ~ T. Such a particle would fall in the first case or in an intermediate case, and may play 
a role in determining the character of the phase transition [52| • However, such a tuning may 



induce unwanted minima in the scalar potential [52J, |53j . We are not going to consider this 



possibility here. Accordingly, the cubic term in the effective potential is 



7 E Imp (as) 



12tt 



where the sum is only over gauge bosons, and the factor 2/3 is due to the fact that only two 
degrees of freedom of the massive vector contribute 3 • 



e. 4 terms. The corrections to A depend logarithmically on T, so the effective value of 
A may be regarded to be constant during the phase transition. Furthermore, these corrections 
are of order /?|/647r 2 , so they can be neglected provided that A > h 2 . For simplicity we will 
assume that this is the case. 

Putting all these terms together we see that, under the above assumptions and approxi- 



mations, the free energy density takes the form displayed in Eqs. (Ifi7|fi8j) . with coefficients 
given by Eqs. (JSJ). 



APPENDIX B: FRICTION COEFFICIENT 

In this appendix we make a derivation of the friction exerted by the hot plasma on the 
bubble walls. For that, we must calculate the departure from equilibrium of the phase space 
population functions, Sf in Eq. (|71jh The friction on the wall has been extensively studied 
in the case of the electroweak phase transition 0, 0, [HI 0, EH EH EH- Our aim here is to 
discuss the general dependence of the friction on the particle content of a theory, so we will 
need to use some approximations in order to keep the description as general as possible. 



1. Fluid approximation 

We begin by considering the contribution of particles with p 3> L~ l (L w = wall width), 
for which the background field varies slowly and the semiclassical (WKB) approximation 
is valid. Since in general L" 1 ^> T, this condition is satisfied for all but the most infrared 
particles [3], which we study below. We follow Refs. 0, E|, but we use a simpler ansatz for 
the deviations from equilibrium distributions. This will suffice for our purposes. We assume 
that the population density of a particle species in the background of the domain wall (that 
moves along the z-direction) is governed by the Boltzmann equation 

[d t + (d p E) d z - (d z E) d Pz ] f = -C [f] , (Bl) 

where E = \J p 2 + m (z, t) 2 is the particle energy, d Pz E = p z /E is the particle velocity, 
—d z E = —d z (m 2 ) /2E is the force on the particle, and C [f] is the collision integral. 



We use the ansantz / = f (E/T — 5), where 



/o (*) = TTT> (B2) 
e x ± 1 

so the deviation from /o {E/T) is 5/ = — /q (E/T) 5. Thus we obtain an equation for 8 by 
linearizing the Boltzmann equation. Keeping only terms of order (m/T) 2 we have 

{^f &tm2 ~ ^ + I 925 ) + ^ [/] = °" (B3) 

The mass of the particle is a function of z — v w t, and so is the perturbation S if we assume a 
stationary state. Thus we make the replacements dtm 2 = —v w (m 2 ) and dtS = —v w 5', where 
the prime means derivative with respect to z. We further simplify Eq. ()B3|1 by making the 
integration J d 3 p/ (27r) 3 . We obtain 

c 2 v w 5' -T5 = dv w m 2 '/2T 2 , (B4) 

where c\ and c 2 are defined by the integrals 

= f d3p f = L f _fp_f 

Cl ~ T*J ( 27 r)V ' ° 2 ~ TiJ (27 r) 3/0 ' 
and we have written the collision integral in the form 

d 3 P C[f] 



(2tt) 3 T 2 



TTS. (B5) 



To lowest order in m/T we have C\f = log2/27r 2 and c^f = 1/12 for fermions, and c^b = 
log (T/m) /2n 2 and c^b = 1/6 for bosons. 

For each particle species i we have a fluid equation of the form ()B4|) . These equations 
are coupled through the collision term (|B5|) . and T is in principle a matrix with indices 
running over all particle species 8 . However, only particles with large Yukawa couplings are 
relevant for the friction force, since they have stronger interactions with the bubble wall. In 
accordance with 0, 0|, in this appendix we call "heavy" these particles with large hi. Notice 
however that heavy particles with large mass fii in the unbroken phase may be thermally 
decoupled and not contribute to the friction at all. The remaining "light" particles can be 
treated as a common background perturbation 5^- The fluid equation for the background is 
simpler and can be solved to eliminate Sb g - Moreover, the heavy particles primarily collide 
with the light particles, so direct coupling between heavy species can be neglected. The effect 
of the light background, though, once eliminated from the equations, is to introduce a weak 
coupling between heavy particles. As a consequence, heavy particles are only weakly coupled 
through the background, and the non-diagonal terms of T are suppressed with respect to 
the diagonal terms by a factor l/g*ii g ht, where g*ii g ht is the number of light species of the 



In Refs. |2ll3| a more complex approximation was made, where the perturbation 5 is split up into three 
different perturbations, S = n/T + E5T/T 2 +p z v/T. In that case there are three fluid equations for each 
particle species, with different rates T M , IV and T„, and there arise additional constants C3 and C4. Our 
approximation corresponds to considering only the term fi/T. 



background (it is proportional to the heat capacity of the plasma). We will therefore neglect 
non-diagonal terms of T in our analysis. Calculating the rates T is well beyond the scope 
of this work. They are of the form a 2 log (1/a) T, where a is a gauge coupling |2j. We will 
assume that in general, T < 10 -1 T. 

The r.h.s of Eq. f|B4j) is the source term of the equation. It is localized at the bubble wall, 
so we expect the same for 5. Therefore we have 5'/ 5 ~ 1/L W . Normally, L w > lOT^ 1 , so if 
the wall velocity is small, the first term on the l.h.s of Eq. (|B4j) is much less than the second 
one and can be neglected. With this approximation the equation has a simple solution, 

S = -v w ^ r . (B6) 
If we now insert 5 fa = —f' Q (Ei/T) Si in Eq. (fTTj) we obtain 

T 2 dm 2 

(vl - 1) <P" + V (0, T) + L- J2 Wii^-St. (B7) 



Replacing the value of 5( given by Eq. ()B6|) . multiplying times 0', then integrating with 
respect to 2, we get 

-V (4> m , T)=v w Y J 9rj - ^ dz. (B8) 

The l.h.s. is the pressure difference between the two phases. It is equilibrated by a friction 
force of the form t]v w . The friction coefficient is thus 

i= u heavy" 

The coefficient c\ for bosons depends on rrii, but it is easy to see that its variation with z 
can be neglected, and we can make the approximation 

ci 6 = logV727r 2 . (BIO) 

To evaluate the integral in ()B9|1 we use the thin wall approximation (|H7j) . 



(.2 1/2, 



dz= / <j) 2 V2Vd(j). 



It is clear that this integral goes like 0^cr. For the model (|BH|). the integral is easily calculated 
using Eq. (|H4*jl . It gives (3/10)0^0", so the friction coefficient is given by 

1 ^ \ 2tt 2 J (Ti/IO^T) T K J 

where Xi = 2 for fermions and Xi = h^ 1 for bosons. According to the arguments above, we 
will make the assumption that the parenthesis in the denominator of the last equation is 
roughly ~ 1. This gives Eqs. (fTHj) and (|02} - 



2. Infrared bosons 



It has been shown 54] that coherent gauge fields can have important contributions to 
the friction. Following Ref. we will estimate the contribution of a gauge boson to 77. 
Infrared boson excitations must be treated classically pq : furthermore, the dynamics of 
the soft fields is overdamped by hard particles jHf|. As a consequence, the equation for the 
population function is given by 

^ d i = -E 2 f + noise, (B12) 
8p at 

which comes from a similar equation for the amplitude of the field. Here, is the Debye 
mass, given by m 2 D ~ gh 2 T 2 , where, according to our previous notation, h is the gauge 
coupling, and g is roughly proportional to the number of particles that couple to the gauge 
field. Averaging over the noise, we get the restoring term —E 2 Sf in the rhs of Eq. (|B12|) . 
Since / = f (E/T) + 5f, and 5f <C / for small v w , we can write 

Inserting in 1)71)1. multiplying by <f)', and integrating as before, we obtain 

t r (rp\ ST V ^ m D [ A ( 2/\2 f d 3 P f' Q 

y ( T ) ^ ~ 2^ 9b^^~ / dz{m h ) / 77-37-™- (B14) 



gauge 



(2tt) j 4pT£ 4 



Since the momentum integral is infrared dominated, we can approximate f' Q (x) ~ — 1/x 2 , 
so the momentum integral yields T/327T 2 mf. With ml = hl<f) 2 , we have 

gauge 

The last integral can be calculated using again the thin wall approximation, Eqs. l)H7)) and 

n dd) r— 2 f 1 , 1 - x . 

There is a logarithmic divergence that must be cut off where the approximations used in 
this derivation break down j54|. Perturbation theory breaks down when mj (0) ~ h 2 T, i.e., 
at (j)/(j) m ~ hbX/E. The kinetic theory description that leads to Eq. ()7T)l breaks down 
when m& (0) ~ L^ 1 , i.e., at d/(j) m ~ \/A//ij. The latter occurs first, so the log is cut off 
at nib ((f>) L w ~ 1- In Ref. j^J] it is argued that the contribution of very infrared degrees 
of freedom is subdominant, since their wavelength cannot resolve the thickness of the wall. 
Hence, the integral in (jB16j) gives to leading log, log0 m /0 = log (mj (0 m ) L„), and the 
friction coefficient is 

gauge w 



which gives Eqs. (fTUj) and 
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